VNIIA-CFAR- 01/12 

Holographic model of 5^ multiband superconductor 

A. Krikun, V. P. Kirilin, and A. V. Sadofyev 
All-Russia research institute of automatics (VNIIA), 
Center of fundamental and applied research, 
Moscow, Russia 
and 

Institute for Theoretical and Experimental Physics (ITER), 

Moscow, Russia 

Abstract 

We construct the holographic model of multiband superconductor. This system is a candidate to 
explain the anomalous features of the iron-based superconductors (e.g. LaFeAsO, BFe2As2, and other 
pnictides and arsenides) . We study the framework which allows the formation of the sign-interchanging 



order parameter. We also calculate the electric AC conductivity, which exhibits a mid-infrared peak, 
related to the interband interaction. 
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I. INTRODUCTION 



The recently discovered family of superconductors, based on the iron compounds {LaFeAsO, 
BaFe2As2 etc.) has attracted lots of attention and revived interest to the unconventional su- 
perconductivity. The question about the possible mechanism of superconductivity is intimately 
related to the symmetry of the order parameter in the system and there is still no consensus 
about the form of this symmetry in iron superconductors These compounds all share 

characteristic multiband feature. Namely they possess multiple conductive bands, which cross 
the Fermi energy level along several Fermi surfaces. In the DFT calculations and ARPES ex- 
periments it was shown, that iron-based superconductors have at least a pair of electron and 
hole Fermi surfaces in the middle (M point) and in the corner (F point) of the Brilluin zone, 
respectively ^. Such a band structure allows the formation of different superconducting conden- 
sates on different conductive bands, what leads to the wide variety of possible symmetries of the 
superconducting order parameter. Among them one of the simplest possibilities was considered 



m 



If the material has two different conductive bands (populated with holes and elec- 
trons, respectively), it is possible that two charged isotropic (S-wave) condensates are formed 
on them, Ai and A2. There are possibilities, where they have the same {sign{Ai) = sign{A2)) 
or opposite {sign{Ai) = —sign{A2)) signs and these arrangements are called S~^^ and order 
parameters, respectively. The two-band system can be studied within the BCS, Eliashberg or 
Ginzburg-Landau approaches [SKZl and one can show, that the 5"^ superconductivity has a 
number of interesting properties. Namely it can be realized even if the intraband interaction is 
weak and the interband interaction is repulsive. In the case of different signs of order parameters 
at different bands the interband interaction may be mediated by the spin fluctuations, and thus 
the mechanism of the Cooper pair creation may be magnetic rather then phononic. Moreover, 
one should anticipate the rich spin dynamics in iron-based superconductors, as they exhibit 
antiferromagnetic instability just near the superconducting region of the phase diagram. The 
strong interband interaction in these systems is also usually related to the anomalous features 



in the AC (optical) electrical conductivity - the peaks in mid-infrared region |9|-|ll|. 



^ There are however exclusions, such as KFeSe, which has only electron-type Fermi surfaces 
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Since the underlying dynamics of the superconducting system under consideration is presum- 
ably strongly coupled, it is reasonable to try to construct a holographic model of the multiband 
superconductor with symmetry of the order parameter, which might describe the variety 
of experimental data and provide some insight into the relations between different features of 
unconventional superconductors. One of the first holographic model of superconductivity has 



been proposed in 



12 



-|l4|. It described the one-band iS-wave superconductor and showed the 



effectiveness of the method by correctly reproducing results of the Ginzburg-Landau theory and 
predicting the nontrivial dependence of AC-conductivity on the frequency. Later other models 
have been constructed for P-wave |l5| and D-wave |l3, 13 superconductors, describing the 



and impurities 



vortex and droplet formation lol and taking into account the effects of crystal lattice |20l. |21| 



22| ^. In our present study we will also use the approach similar to 28|, where 



the multiband superconductor was considered, although our construction is different from that 



used in 



28 



providing us with the possibility to treat the phases of condensates and electric 
current straightforwardly. The paper is organized as follows: in Section II we introduce the 
holographic model of superconductor, in Section III we study the phase diagram of the 
model and the features of different condensate patterns. Section IV is devoted to the calcula- 
tion of the electric conductivity of the dual superconductor which possesses anomalous features 
related to the interband interaction, the conclusion is given in Section IV. Appendix A is de- 
voted to the derivation of the full set of equations of motion in our model. In Appendix B we 
elaborate the numerical equation of state. 



II. HOLOGRAPHIC MODEL 

We start building the holographic model of 5*^ superconductor by identifying the degrees 
of freedom in the bulk theory. First of all, the S symmetry of the order parameter suggests. 



that its dual bulk field should be a scalar 



12 



13| . Moreover, in the two-band material, we 



are considering, there are two possible gaps Aj(i = 1, 2) in different bands, as we pointed out 
earlier, each with its own complex phase. This allows us to argue, that there may be the 



^ Reviews of the subject include |23l427| 
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hidden global symmetry between bands in the boundary theory. The use of such a symmetry 
is the main ingredient of the our model, which leads to the results, qualitatively similar to the 
experimentally observed phenomena. Thereby we organize the two gaps in the fundamental 
representation of U{2) group. This global symmetry in the boundary theory leads, according 
to the general holographic principle to the similar gauge symmetry in the gravity side of the 
correspondence. Thus, in the bulk we obtain the fundamental scalar field = (p2)'^, charged 
under the nonabelian U (2) gauge symmetry. The method of describing multiband holographic 



systems by use of the hidden global symmetry was adopted also in [28|, but the usage of U{2) 
group instead of SO (3) allows us to define phases of complex condensates separately, what is 
otherwise rather tricky. The action of the bulk theory is 

1 



S = d xdr 



(1) 



where F^^, = d^Ay — d^A^ — i[A^,Au], A^ is the U{2) gauge field, the covariant derivative is 
-D^0* = (9^0* + i{Afj,(j)Y and we rescaled the bulk gauge coupling to unity. We choose the mass 
of the scalar field according to the single band case [3]: m| = —j2- We will use the probe 
approximation, neglecting the backreaction of the bulk fields on the gravity, so the background 
metric is just the AdS-Schwarzshield planar black hole: 

ds'' = -f{r)dt^ + j^ + r\dx^ + dy^), (2) 

2 

The blackening factor /(r) = ^ — ^ defines the position of the black hole horizon Vh, which 
is related to the temperature of the boundary theory as 

In the paper we will always rescale the length units so that L = 1. 

In what follows we try to find the solution to equations of motion derived from ([T]), which 
describes the condensation of the charge carriers in both bands of the boundary theory. In 
accordance to the holographic dictionary the condensate A of the operator is proportional to 
the normalizable mode in the r — )> oo expansion of the classical solution for the corresponding 
bulk field (p 

A(l) (/,(2) 

<P{r) =^ + ^ + ... (4) 



In the present case both modes are normahzable, and we choose the second one 0^^^ = A, thus 
defining a particular boundary system to be described (see 12|, |l3| for discussion of possible 



choices). The first mode (f)^^^ is related to the source of the operator in the boundary theory 
and should be put to zero, as there are no such sources in the problem under consideration. To 
find the solution with two condensates we use the ansatz: 

f cos(9)\ fl o\ [o l\ 

0(r) = Mr) \, A = M{r) \ dt + A(r) dt, (5) 

Vsin(^)/ V° V VV 

with constant 9, which describes various possible relations between two condensate. We find, 
that to satisfy the equations of motion of the model, in given ansatz 6 must assume values 

0+ =(4n + 1)^, (n G Z), S++ state (6) 

e- =(4n + 3)^, G Z), state (7) 

Actually the state which the scalar field chooses to condense in depends stronlgy on the nonan- 
belian component of the gauge field, turned on in the particular ansatz, as it is proportional 
to corresponding eigenvector. We will comment on other possible solutions bellow. Here the 
6 phases with superindex + and — describe the S^^ (A, A) and S"^ (A, —A) symmetries of 
the order parameter, respectively. The time component of the abelian gauge field M(r) is dual 
to the number of particles operator in the boundary theory c|cj, hence its boundary value is 
equal to the corresponding chemical potential. In this ansatz we consider the equal chemical 
potentials on two bands (//i = ^2 = A*), assuming, that the electron and hole bands of the 
superconductor are doped by the same amount of carriers. It is straightforward to generalize 
this model to the case of different chemical potentials in the bands by introducing the third 
component of the nonabelian gauge field {A^ ~ diag{l^—l)), but we won't do this here for 
simplicity. The off-diagonal component of the nonabelian gauge field A is also dual to the 
scalar current in the boundary theory, similar to the particle number operator, but this time 
it describes the jumps of particles from one band to another: c\c2 + c\ci. The presence of the 
source for such an operator would mean that there exists an external potential which results 
in the mixture of different bands, so that only mixed states which diagonalize the matrix of 
"chemical potentials" , have definite particle numbers. 
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In Appendix A we check that the ansatz satisfies the full set of equations of motion 
following from ([T]) if the functions <po{r), M(r), A(r) are solutions to (primes denote derivatives 
with respect to r) 



M 
A 



f 



(M ± A)^ 

7^ 



fir) 



(8) 



A" 



r 



-A't^(M±A) 
r f 







Here and below upper signs correspond to S~^~^ ansatz and lower - to S^. One can combine 
the equations for M and A in such a way, that the system will look exactly like the equations 
of motion for a single scalar, studied in [12| 

(M±A)2 2 



^0 • 



M±A 



(M±A)" + -(M±Ay 
r 



0, 



(9) 



2ji 
f 



{M ± A) = 0. 



This system describes a scalar field, that gets a nonvanishing normalizable component at low 
temperatures T < T^.. Composing the another linear combination of equations dH]) one gets 



MtA: (MtA)" + -(MtA)' = 0. 



(10) 



In 



12| the critical temperature Tc was related to the density of charge carriers, or the 



normalizable component of the field (M(r) ± A(r)). In our case the charge carriers density p 
and normalizable component of (M(r) ± A(r)) are generally not equal, so we should study the 
equation of state for our system in more detail. First of all let us fix the notation, denoting the 
asymptotic coefficients for the fields as 



A(r) 
M(r) 



r 
P 

/i- -. 

r 



(11) 



Here p is proportional to the charge carriers density (cjcj) as by holographic prescription the 
normalizable mode of the field corresponds to the mean value of dual operator. Similarly, J is 
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FIG. 1: The current J = (cjc2 + CgCi) in S~^~^ condensed state of two-band superconductor 
(in S"^ it has opposite sign) in absence of the corresponding source A. 

proportional to the value (c|c2 + ^ci) which can be interpreted as the particle jumping rate 
between the bands. 

As the system of equations (p|). (fTO!) is scale invariant, we can choose one of the parameters 
of our model to define the scale of all dimesionful quantities. Taking = /iq ~ const would 
mean we are considering grand canonical ansamble and the number of particles can change with 
temperature. The second possibility is to take p = po ~ const as reference scale, considering 
the canonical ansamble with constant number of particles on the bands. In what follows we 
assume the latter point of view. Moreover, we find it convenient to keep constant the ratio of 
/i and A in order to control the mixing angle of band states. 



Equations fl9|). f|T0|) unambiguously define the state of the system. The procedure of getting 
the corresponding parameters from the numerical solution is described in Appendix B. In the 
end of the day we find, that the solution with nonzero scalar condensate exists in both ansatze 
|5] below certain critical temperature, which is ( 1371) 



A 

a = — ~ const 



(12) 




for S^^ state 



(13) 



Tc ^ 0.118v/po(l - a) for state 
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It is also noteworthy, that in the condensed phase of our two-band superconductor there is 
a nonzero current J = {c\c2 + c^ci) even in the absence of external source A = 0, a = (see 
Fig. [1]). This can be interpreted by noting that strong interband interaction is mediated in the 
bulk and leads to the interband jumping of the charge carriers. 

III. PHASE DIAGRAM 

One can see that in the situation when the off diagonal component of chemical potential 
A is absent, the two possible configurations of the condensates are degenerate: they have the 
same critical temperature fll3p and the same free energy, as the equations for different 
ansatze switch places when the sign of A is changed. Consequently, any perturbation of the 
system would brake this degeneracy and we can consider nonzero A as one of the possibilities. 
One should note, however, that if A = then there is no reason the ansatz ([5]) would have 
minimal free energy among other possible solutions. Indeed, in such a situation the solution 
with nontrivial bulk profile of any nonabelian gauge field component should have the same 
free energy as ([5]) by the symmetry reasons. But as we have already noted, the condensation 
pattern of the scalar field (and consequently the symmetry of the superconducting gap) is 
crucially dependent on the type of condensed gauge component, as it is proportional to one 
of its eigenvectors. Namely, in order to exhibit the gap in both conductive bands the system 
should contain a perturbation, which leads to the condensation of component of the gauge 
field and turning on the potential A is a reasonable choice to achieve this effect. In this section 
we will check the free energy of different condensation patterns at nonzero A to find, which 
solution is the stable one. 

First of all we note, that the critical temperatures of S~^~^ and differ at a 7^ 0. As we 
would like to study the condensation, we will take in what follows a = —0.05, corresponding 
to small negative A. Hence on the plot of the order parameter value with respect to temper- 
ature (normalized to the critical temperature Tq at A = 0) Fig|2](a) we see, that the critical 
temperature of S'^'^-phase in this situation is higher. 

To find out, what configuration of the order parameter is stable we compute the free energy 
F of our solutions. By the holographic principle, F is related to the on-shell action for given 
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FIG. 2: Dependence of the condensate 
value on the temperature for S"^ state 
(purple curve) and S""^"*" state (brown 
curve) at negative A = — 0.05/i. Black 
dashed curve - degenerate state at A = 0. 




T 

1.4 r„ 



FIG. 3: Free energy of S"^ state (purple 
curve) and S^~^ state (brown curve) at 
A = —0.05/i. Blue curve - free energy of 
the normal (nonsuperconducting) phase 
at the same A. 



solution. The action on our ansatz 



IS 



on~shell 



d xdrr 



{drMf + {drKy 



/(5r0o)' + y0o(^±A)' + 202 



Integrating by parts and using the equations of motion we can rewrite it in the form 

„2 



on— shell 



V2 / dt 



[MdrM + AdrA 



odr(po} 



rh 



00 iM±Ay 



fir) 



here V2 is 2d volume. This yields for the F = TS, 



on— shell 



Vo 



{fip + A J) 



<Pl (M ± A)^ j^^dr 



(14) 



(15) 



(16) 



We can clearly interpret the first term as the value of F in the normal, uncondensed phase. 
The second one is always negative, rendering the superconducting phase thermodynamically 
favorable. We also observe that the condensate with higher for given A is the one which 
persists when we further lower the temperature, since it has lower F. That is, for A > 0, 
T^^ > and F^^ < F^, and vice versa for A < 0. Indeed we see the corresponding behavior 
on the numerical plot of F of normal and condensed phases on Fig|3l 
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IV. ELECTRICAL CONDUCTIVITY 



To compute the electrical conductivity in the condensed phase of holographic supercon- 
ductor we first of all identify the elecromagnetic current in the model. As the two conductive 
bands are populated with the carriers with opposite charges (electrons and holes), the electro- 
magnetic field will couple differently with each of the condensates. Namely, the electromagnetic 
current in this model is dual to the nonabelian component of the gauge field 




/13(1) 

^3(0) ^ ^ (17) 



Thus to study the electric conductivity we consider the perturbation of ^4^ on the background of 
the condensed solution and calculate the relation between its normalizable and non-normalizable 
components, the former related to the electric current in the boundary theory and the latter 
to the vector potential of external electric field: 



Let us turn on the electric field along the x-axis. By looking at the equations of motion (see 
Appendix) one can notice, that multiplied by A] = A(r) provides a linear source for the 
other gauge field component A^. In its turn enters the equation of motion for A^. Thus 
about the background ([5]) we get the system of coupled linear equations on A^jA"^ ~ e'"^*, 
describing the perturbation of the electric field 

Al : d'rAl + ^drAl + - 7 + + = 0' (19) 

Al : d^.Al + fjd^Al + - I + 4^) - ^^Al = 0. (20) 
Which can be readily rewritten as a couple of equations on the complex functions A"^ = Al±iA'^ 



+ Q^drA^ + 
fir) 



f{rf f{r) 



^± = 0. (21) 



Now to compute the conductivity we solve these equations with outgoing wave boundary condi- 
tions on the horizon: A'^,A^ ~ (r — r/i)^'h , and equal sources on the AdS boundary (r — )■ oo): 
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T 



(a) (b) 

FIG. 4: The real (a) and imaginary (b) parts of AC conductivity depending on frequency at 
negative A = — 0.05/i: purple curve - S*^ state, brown curve - S''^+(this state is quasistable at 
given A). Blue curve - degenerate case with A = 0. Dashed black curve - the ordinary 
one-band superconductor with the same parameters. Temperature is constant T ^ O.STq 

= A^^^\ Then we compute the relation of the normalizable and non-normalizable 
modes of Al{r) = i(^+(r) + ^~(r)). 

The dependence of the AC conductivity on the frequency for constant temperature 
(T ^ O.8T0) for the S± and S^~^ ansatze is plotted on FigJH The imaginary part of con- 
ductivity exhibits a pole at a; = in all cases, what is a signal of superconductivity and is 
related to the 5-function at a; = in the real part of a by the Kramers- Kronig relation. This 
pole ensures us, that our model indeed describe a superconductor. More interestingly, even 
without introducing the external interband potential A = (blue curve) we see drastic change 
of the behavior of the real part of the conductivity in comparing with the ordinary one-band 
superconductor with the same parameters (dashed curve). When we consider finite off-diagonal 
chemical potential A = — 0.05/x the peak in the conductivity lowers for S"^ superconductor (pur- 
ple curve) and rises in the 5'"'""'' superconductor (brown curve). One should note however, that 
S~^~^ state is quasistable at such A so in the real material the state with lower peak will survive. 
We can conclude, that the mid-infrared peak in the AC conductivity should be observed in 
any multiband system, regardless of the relative sign of the order parameters, provided the 
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off-diagonal chemical potential (which can be related to the density of impurities, see below) is 
not too large. At large A the peak disappears, but the hump at ^ ~ 18 becomes larger. 

The peak in the mid-infrared region of frequency is usually related to the interband scattering 
and was obtained, for instance, in the calculation of the conductivity of state within the 
Eliashberg theory framework Moreover, in {g] it was discussed, that the interband scattering 
can be controlled by the impurities and this statement is in strong relation with our model 
in case of nonzero A. The effects of impurities can be described in holographic models by 
introducing external sources for vector operators 22[, so we find that introducing nonzero A 
can be treated indeed as an impurity effect. 

One can also find such peaks in the experiments on measuring the optical reflectivity on the 



samples of LaFeAsO 



le expei 



10|, BaFe2As2U\ 

del shows, that the mid-infrared 
peak in AC conductivity is a consequence of the multiband structure of superconductor and 
doesn't need any additional mechanism to be explained. 

In the end of this section we should underline the importance of the different sign of electric 
charge of carriers (electrons and holes) at different conductive bands of our superconductor. 
Indeed if the carriers had the same charge sign on both bands, the electromagnetic current 
would couple to the Abelian part of the gauge field and would not interact with other nonabelian 
components in the bulk. In such situation we would obtain an equation for electromagnetic 
fluctuation, equivalent to the case of ordinary one-band superconductor |l2| and would not see 
any interesting features in the AC conductivity (it would look like the dashed curve on Figj4]). 
Thus the fact that one band is hole-type and another is electron-type is crucial. 



V. CONCLUSION 



In this work we constructed a holographic model of two-band isotropic superconductor. 
The model includes a charged scalar field in the fundamental representation of U{2) gauge 
group, whose vacuum expectation value describes the condensates on the different bands of 
the superconductor. We found that the pattern of condensation is crucially dependent on 
the perturbation, that resolves the degeneracy of all possible ways of spontaneous breaking 
of the gauge symmetry by scalar field. We observe the condensation with equal gaps on two 
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bands if the effect of perturbation can be expressed as nonzero off-diagonal chemical potential, 
describing the interband interaction. Moreover, we find that the final state, in which the 
system condenses (with the maximal value of the critical temperature) depends on the sign of 
this interband interaction and if the off-diagonal chemical potential A is negative, the state 
is more favorable. This allows us to argue that the next step in development of holographic 
model would be to introduce some additional ingredients (possibly the magnetic subsystem, 
described by the spin fluctuations or impurities) in order to generate such an interaction (off- 
diagonal chemical potential) and resolve the degeneracy without the need in external potentials. 

In our study of the AC electric conductivity we find, that the mid-infrared peak in the real 
part of the spectra is a consequence of the multiband nature of the material and is saturated 
by the interband interactions if the carriers on different bands have opposite electric charges 
(electrons and holes). As our analysis is performed in the probe approximation (the backre- 
action of the metric is not taken into account) we can not study the dependence of various 
features of this peak on the temperature at large enough temperature range. The account for 
the backreaction would also help to move to lower temperatures, where the gap formation can 
be seen in the conductivity spectrum. The study of these features would shed more light on 
the nature of interband interaction in superconductors, so we leave this subject for future 
work. 

Our model gives results in close relation with the similar studies within the Eliashberg 
framework and it would be interesting also to explore nontrivial features found there in this 
model, for instance the transition from order to S~^~^ order when the impurity scattering is 
strong. 
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VI. APPENDIX A: EQUATIONS OF MOTION 

In appendix we consider the full set of equations of motion of our model ([1]) and check, that 
the ansatz (j5j) satisfies them. First of all introduce the parametrization of the scalar field in 
the fundamental representation of U{2) 



Ooe 



sin(6')e*'^ 
cos(6')e~*'^ 



and the U{2) gauge field 



— ^fj. 




The action of the model ([T]) in this notation takes the form 



S 



d xdr 



2^" 



- {d,<Po? - <Pl{A, + d,aY - 0^(4 + a,7)' 
+ 2<Pl cos{2e)g^'' (4 + d^a) {Al + d,-f) 

- 2<Pl sm{2e)g^'' (i^ + d^a) {Al cos(27) + Al sin(27)) 

- 20^(4 sin(27) - Ajcos(27))a^^ - 4>l{{Alf + {A 



\2\2 



(22) 



(23) 



(24) 

(25) 
(26) 
(27) 
(28) 
(29) 

(30) 



where F = dA and = dA\ From this action we derive the equations of motion in arbitrary 
gauge 
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--d„ 



'-9 



'-9 



'-9 



- 2 cos{2e)g^'' (i^ + d^a) {AI + d,^) 
+ 2 sin(2^)(7'^'^(4 + d^a) [Al cos(27) + sin(27)) 
+ 2(4 sin(27) - Al cos{2^))d^e + ((4)^ + {A 



If 



+ (Pi sm{2e) {Al cos(27) + Al sm{2^)) 



--da\/^9 



(j)l{Al + a,7) - (Pi cos(2^) {A, + d,a) 

2(Pl sm{2e)g^'' (i^ + d^a) {Al sm{2^) - Al 005(27)) 
2(Py{Alcos{2^) + Alsm{2^))d,e = 



99 



(PldJ + (Pl{Alsm{2j) - Alcosi2j)) 

- 2(Pl ^in{2e)g^'' (i^ + d^a) {AI + ^,7) 

- 2(Pl cos(2^) (7^'^ (4 + a^a) (Aj, 005(27) + Al sin(27)) = 

~99'''9''^Fxp-cPl9''{K + dp^) 

+ (t>lcos{2e)g^^{Al + d,^) 

- 02sin(2^)^^''(y4jcos(27) + Ajsin(27)) = 

K + 2(^a4 - ^IK)] - ^9''9''{F^,Al - F^^Al) 



99''9''' 



Ag'^Pg 



(Pi sin(2^) cosi2^)g>^P{A, + d,a) - (Pi sin (27) (/'^"^.^ - (P^^Al = 
Fl + 2(A^,A;, - AlAl)] - 2g^''g^\Fl^Al - F^^Al) 



^9^'9 



- Al{A\Al) 



-99''9''^ 



(g sin(2^) ^m{2^)g^^{Ap + dpo) + (g cos{2^) g^^ d ^9 - (Plg^'Al = 
Fl + 2(A1AJ - A^4)j - 2^/''''^^^(F^2,Ai - F^Al) 
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+ (fl cos{2e)g^^P{Ap + a,a) - (gg^'^Al + ^,7) = 



One can check, that using the ansatz the equations of motion reduce to 



a 



A, 



A]. 



.d,V^g^''<l>l{A,±Al)=0 
--duV^g^'g^'F^p - (f^lg^'^^A, ± Aj) = 
--d^V^g^'^g'-'Fl T 49''{a, ± 4) = 



And equations for 7, 9, A^, A^ are zero. The equation of motion for a follows from equation for 
A, what is an evidence of the gauge symmetry. We see, that our ansatz ([5]) satisfies the general 
equations of motion if the functions (f)o{r), M{r) = At{r),A{r) = Al{r) are solutions to ([8]) 



VII. APPENDIX B: NUMERICAL EQUATION OF STATE 



Solving the system ([9j) numerically we get a one parameter family of solutions 

{B''{r),(f)\r),rl}, where 



BUr) = MUr)±AUr) 



J \2 



(j)\r) > 0* 



(31) 
(32) 



On the other hand one can easily find the solution of ffTOl) satisfying the boundary conditions 
(ITT]) and vanishing at the horizon (as due to the divergence of g^* at r = we need to impose 
the condition Atijh) = 0) 



1 _ ^ 

r 



(33) 



These relations allow us to derive the corresponding parameters of the z-th solution. In the 
S~^~^ case they are 



2po 



1 



A* 



/i 



2po 



(34) 
(35) 
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The case 5'''= differs by the change of sign in front of a and J. The critical parameters of the 
numerical solution are 

b^"- ^ 6f ^ 4.06 (36) 
Thus we get the critical temperatures for S~^~^ and S"^ states are 



2po 



47rY 4.06 + 4.06 (i^) 



2po 



47r y 4.06 + 4.06 {^) 



0.118Vpo(l + a) 



0.118v/po(l - a) 



(37) 



(compare with the result of [l2|: T ^ 0.118^/po 
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